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This paper presents a complete description of noise-induced decay of a metastable state in a 
wide range of noise intensity. Recurrent formulas of exact moments of decay time valid for arbitrary 
noise intensity have been obtained. The nondecay probability of a metastable state was found to 
be really close to the exponent even for the case when the potential barrier height is comparable or 
smaller than the noise intensity. 
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I. INTRODUCTION 

Overdamped Brownian motion in a field of force (Markov process) is a model widely used for description of noise- 
induced transitions in different polystable systems. In many practical tasks (e.g. tasks of Josephson electronics, 
kinetics of chemical reactions and so on) it is enough to know the probability of transition or even only a time scale 
of transition. When the transition occurs over a potential barrier high enough in comparison with noise intensity, 
the probability of transition is a simple exponent ~ exp(— t/r) [JjJ, where r is the mean transition time. In this case 
the mean transition time gives complete information about the probability evolution. The boundaries of validity of 
exponential approximation of the probability were previously studied in ||. In authors extended the Mean 
First Passage Time to the case of " radiation" boundary condition and for two barrierless examples demonstrated good 
coincidence between exponential approximation and numerically obtained probability. In a more general case the 
exponential behavior of observables was demonstrated in || for relaxation processes in systems having steady states. 
Using the approach of " generalized moment approximation" the authors of || obtained the exact mean relaxation time 
to steady state and for particular example of a rectangular potential well demonstrated good coincidence of exponential 
approximation with numerically obtained observables. The considered in || example of the rectangular well does not 
have a potential barrier, and the authors of that paper supposed that their approach (and the corresponding formulas) 
should also give good approximation in tasks with diffusive barrier crossing for a wide range of noise intensity. 

In the frame of this paper we consider a different case than in ||, S: we present investigation of nondecay 
probability of a metastable state. We treat the decay as a transition of Markov process trajectory outside the region 
of a metastable state. We consider namely the probability to find a relization of Markov process in a given interval, 
but not the probability to pass the boundary for the first time (First Passage Time formalizm), since usually in 
experiment we can measure only the first one. In most of applied tasks no absorbing boundaries can be introduced 
and the precision of the used equipment allows to register the only fact of leaving the considered domain, but the 
boundary may be crossed many times during the transition (infinite number of times for Markov process, since this 
process is not differentiable) . 

Using the approach proposed by Malakhov Q , Q , that requires only knowledge of the behavior of a potential at 
±oo, we have decomposed the nondecay probability into a set of moments (cumulants), obtained recurrent formulas for 
these moments and approximately summarized them into the required probability. The obtained nondecay probability 
demonstrates exponential behavior with a good precision even in the case of a small potential barrier in comparison 
with noise intensity. 

II. MAIN EQUATIONS AND SET UP OF THE PROBLEM 

Consider a process of Brownian diffusion in a potential profile <P(x). Let a coordinate x(t) of the Brownian particle 
described by the probability density W(x,t) at initial instant of time has a fixed value x(0) = xo within the interval 
(c, d), i.e. the initial probability density is the delta function: W(x, 0) = S(x — xq), xq € (c,d). 
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In this case the one-dimensional probability density W(x,t) is the transition probability density from the point Xq 
to the point x: W(x,t) — W{x, t; xq, 0). It is known that the probability density W(x, t) of the Brownian particle in 
the overdamped limit satisfies to the Fokker-Planck equation (FPE) : 



dW(x,t) _ 8G(x,t) 1(8 



dt dx B I dx 
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with the delta-shaped initial distribution. Here B = h/kT, G(x,t) is the probability current, h is the viscosity (in 
computer simulations we put h = 1), T is the temperature, k is the Boltzmann constant and f(x) = <P(x)/kT is 
the dimensionless potential profile. In this paper we restrict ourselves by the case of metastable potentials, i.e. we 
consider an overdamped Brownian motion in a potential field (p(x) in systems, having metastable states, such that 
ip{— oo) = +oo and </?(+oo) = — oo. This leads to the following boundary conditions: G{— oo, t) — W(+oo, t) = 0. 
Note, that the results obtained may be generalized for potentials of arbitrary types, e.g. for such that <^(±oo) = oo. 

It is necessary to find the probability P(xo,t) of a Brownian particle, located at the point xo (t = 0) within the 

d 

interval (c, d) to be at the time t > inside the considered interval: P(xq, t) — J W(x, t)dx. Further we for simplicity 

c 

will call the probability P(xo,t) as nondecay probability. We suppose, that c and d are arbitrary chosen points of an 
arbitrary potential profile <p(x) and boundary conditions at these points may be arbitrary: W(c, t) > 0, W(d, t) > 0. 
In this case there is the possibility for a Brownian particle to come back in the interval (c, d) after crossing boundary 
points. 

III. MOMENTS OF DECAY TIME 

Consider the nondecay probability P(xo,t). We can decompose this probability to the set of moments. On the 
other hand, if we know all moments, we can in some cases construct a probability as the set of moments. Thus, 
analogically to moments of the First Passage Time (|]- Q we can introduce moments of decay time T n (c,xo,d) (or, 
generally, moments of transition time, see S, where it was performed for the probability Q(xo, t) = 1 — P(xo, t)): 

oo 

J t ndp(goM dt 

r n (c, xo, d) =< t n >= — ° r — -. (2) 

Here we can formally denote the derivative of the probability divided by the normalization factor as w(xo,t) and 
thus introduce the probability density of decay time w(xo,t) in the following way ||: 

dP(x ,t) i 

WM= dt lP(x Q ,oo)-P(x ,0)Y (3) 

It is important to mention that the moments of decay (transition) time (||) is a generalization of the well-known 
moments of the First Passage Time for the case of arbitrary boundary conditions (see discussion in 0]). For example, 
in the considered case of the potential tp(x) (such that (p(—oo) = +oo and ^(+00) = —00) the moments of decay time 
coincide with the corresponding moments of the First Passage Time, if a reflecting boundary at the point c and an 
absorbing boundary at the point d are introduced. On the other hand, if we consider the decay of metastable state 
as transition over a barrier top, and compare mean decay time obtained via approach discussed in the present paper 
(case of a smooth potential without absorbing boundary) and the mean First Passage Time (MFPT) of the absorbing 
boundary located at the barrier top, we get two times difference between these time characteristics even in the case of 
a high potential barrier in comparison with the noise intensity. This is due to the fact, that the MFPT does not take 
into account the backward probability current and therefore is sensitive to the location of an absorbing boundary. For 
the considered situation, if we will move the boundary point down from the barrier top, the MFPT will increase up to 
two times and tend to reach value of the corresponding mean decay time, which is less sensitive to the location of the 
boundary point over a barrier top. Such weak dependence of the mean decay time from the location of the boundary 
point at the barrier top or further is intuitively obvious: much more time should be spent to reach the barrier top 
(activated escape) than to move down from the barrier top (dynamic motion). 

The required moments of decay time may be obtained via the approach proposed by Malakhov ||, ||. This 
approach is based on the Laplace transformation method of the FPE ([j]) . Following this approach, one can introduce 

the function H(x,s) = sG(x,s), where G(x,s) — J G(x,t)e~ st dt is the Laplace transformation of the probability 



current, and expand it in the power series in s: 
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H(x, s) = sG(x, s) = H (x) + sifi(x) + s 2 H 2 (x) + ... (4) 
It is possible to find the differential equations for H n (x) (see Q, ||; dHo(x) / dx = 0): 

dffi(x) 

— — — = d(x - x ), 

d?H n {x) d<p(x) dH n (x) RTT ... (5) 

— "r~5 1- — j -, =BH n - l (x), n = 2,3,4,... 

Using the boundary conditions W(+oo,t) = and G{— oo,t) = 0, one can obtain from ([|) H\{x) = l(x — xq) and 

H 2 (x) =-B J e-vW Je^yh(y-x )dydv, 

x oo \ u / 

H n (x) = -B J e-vW fevMH^Mdydo, n = 3,4,5,... 

— oo ^ 

Why did we calculate this recurrent formula for the functions H n (x)7 The matter is, that from formula (Q) (taking 
the integral by parts and Laplace transforming it using the property P(xo, 0) — sP(xq,s) = G(d, s) — G(c, s) together 
with the expansion (^)) one can get the following expressions for moments of decay time: 

n (c,x ,d) = -(H 2 (d) - H 2 (c)), 

T 2 (c,x ,d) = 2(H 3 {d)- H 3 (c)), ( . 

T 3 (c,x ,d) = -2 • 3(H 4 (d) - H 4 (c)), ... {l> 
T n (c,x ,d) = (-l) n n\(H n+1 (d) - H n+1 (c)). 

One can represent the n-th moment in the following form: 

T n (c,x ,d) = n!r"(c,xo,rf) + r„(c, x , d). (8) 

This is a natural representation of r n (c, Xo,d) due to the structure of recurrent formulas (|^), which is seen from the 
particular form of the first and the second moments for the case c = — oo (c < xq < d). From the recurrent formulas 
@i (0) one can obtain: 

d oo d x 

Ti (-oo, a*, d) = B I f e-^ x) dx- [ e v{v) dv- [ e"^ 1 ' / e v{v) dvdx\. (9) 



Xq Xq Xq 



j(— oo, x , d) = 2£? 2 |[n(— oo, x , d")] 2 + 



d oo v 

e-*W./e*>/e-*>/e*W„*- (10) 

-oo ceo d 

d a; i) 



Using the approach, applied in the paper by Shenoy and Agarwal [flQf| for analysis of moments of the First Passage 
Time, it can be demonstrated, that in the limit of a high barrier Atp 3> 1 (Aip = A<P/kT is the dimensionless barrier 
height) the remainders r n (c, xq, d) in formula (^) may be neglected. For A<^ m 1, however, a rigorous analysis should 
be performed for estimation of r n (c, Xo, d). Let us suppose, that the remainders r n (c, Xo, c?) may be neglected in wide 
range of parameters and further we will check numerically when our assumption is valid. 

The cumulants of decay time ee n [O] , @] are much more useful for our purpose to construct the probability 
P(xq, t), that is the integral transformation of the introduced probability density of decay time w(xq, t) (||). Unlike 
the representation via moments, the Fourier transformation of the probability density (H) - the characteristic function 
- decomposed into the set of cumulants may be inversely transformed into the probability density. 

Analogically to representation for moments (|J), similar representation can be obtained for cumulants ae n : 

ae n (c, xo,d) = (n- l)!ae"(c, x Q , d) + R n (c,x ,d). (11) 
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It is known that the characteristic function Q(xq,uj) = J w(xo,t)e : ' ult dt (j = y/—l) can be represented as the set of 

o 

cumulants (w(x , t) = for t < 0): 



Q(xq, u>) = exp 



E 

.n=l 



ee n (c,x ,d) 



(12) 



In the case, when the remainders R n (c, xq, d) in ( pi] ) (or r n (c, xq, d) in (JsJ) ) may be neglected, the set ( |l2| ) may be 
summarized and inverse Fourier transformed: 



-t/T 



w(x ,t) 



where r is the mean decay time 0, || (r(c, jgq, d) = t\ = aei): 



(13) 



T(c,x ,d) = B{ / e v(a:) / e~ lp{v) dvdx + / e v(a:) dx / e-^du 



(14) 



This expression is a direct transformation of formula (0), where c is arbitrary, such that c < xq < d. 

Probably, similar procedure was previously used (see~]7) , [jl0| , Q , |l3| ) for summation of the set of moments of the 
First Passage Time, when exponential distribution of the First Passage Time probability density was demonstrated 
for the case of a high potential barrier in comparison with noise intensity. 



IV. NONDECAY PROBABILITY EVOLUTION 



Integrating probability density (|T^), taking into account definition (||), we get the following expression for the 
nondecay probability P(xo,t) (P(xq,0) = 1, P(xo,oo) — 0): 

P(xo,t) =exp(-t/r), (15) 

where mean decay time r is expressed by jl4|) . Probability (|l^) represents a well-known exponential decay of a 
metastable state with a high potential barrierjpj . Where is the boundary of validity of formula ([If]) and when can we 
neglect by reminders r n and R n in formulas (|q),(fi"l|)? To answer this question we have considered three examples of 

d 

potentials having metastable states and compared numerically obtained nondecay probability P(x$, t) = J W(x, t)dx 

__ c 

with its exponential approximation ( |T5| ) . We used the usual explicit difference scheme to solve the FPE (||) , supposing 
the reflecting boundary condition G(cb,t) = (q, < c) far above the potential minimum and the absorbing one 
W(df,,t) = (db > d) far below the potential maximum, instead of boundary conditions at ±oo, such that the 
influence of phantom boundaries at Cb and db on the process of diffusion was negligible. 

The first considered system is described by the potential 4>(x) — ax 2 — bx 3 . We have taken the following particular 
parameters: a = 2, b = 1 that leads to the barrier height A«£ ps 1.2, c = —2, d = 2a/3b, and kT = 0.5; 1; 3. 
The corresponding curves of the numerically simulated probability and its exponential approximation are presented 
in Fig.l. In the worse case when kT = 1 the maximal difference between the corresponding curves is 3.2%. For 
comparison, there is also presented a curve of exponential approximation with the mean First Passage Time (MFPT) 
of the point d for kT = 1 (dashed line). One can see, that in the latter case the error is significantly larger. 

The second considered system is described by the potential <P{x) = ax 4 ~bx 5 . We have taken the following particular 
parameters: a — 1, b = 0.5 that leads to the barrier height A<£ w 1.3, c = —1.5, d = 4a/56, and kT — 0.5; 1; 3. The 
corresponding curves of the numerically simulated probability and its exponential approximation are presented in 
Fig. 2. In the worse case (kT = 1) the maximal difference between the corresponding curves is 3.4%. 

The third considered system is described by the potential •P(x) = 1 — cos(x) — ax. This potential is multistable. We 
have considered it in the interval [—10, 10], taking into account three neighboring minima. We have taken a = 0.85 
that leads to the barrier height w 0.1, c = — tt — arcsin(a), d — tt — arcsin(a), xq — arcsin(a), and kT = 0.1; 0.3; 1. 
The corresponding curves of the numerically simulated probability and its exponential approximation are presented 
in Fig. 3. In difference with two previous examples, this potential was considered in essentially longer interval and 
with smaller barrier. The difference between curves of the numerically simulated probability and its exponential 
approximation is larger. Nevertheless, the qualitative coincidence is good enough. 
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Finally, we have considered an example of metastable state without potential barrier: <P(x) — —bx 3 , where 6=1, 
xq = —1, d = 0, c = —3 and kT = 0.1; 1; 5. By dashed curve an exponential approximation with the MFPT of the 
point d for kT — 1 is presented. It is seen, that even for such example the exponential approximation (with the mean 
decay time (|l4|)) gives an adequate description of the probability evolution and that this approximation works better 
for larger noise intensity. 

V. CONCLUSION 

In the present paper the decay of metastable states, described by the model of Markov process, has been considered. 
Recurrent formulas of exact moments of decay time, valid for arbitrary noise intensity, have been obtained. Some 
concrete examples of metastable states have been analysed numerically, and the time evolution of the nondecay 
probability of a metastable state is found to be really close to the exponent even for the case when the potential 
barrier height is comparable or smaller than the noise intensity if the exact mean decay time ( |l4| ) is substituted into 
the factor of exponent. 

For all investigated examples, the exponential approximation gives an adequate behavior of the probability. This 
approximation may be used in a wide range of parameters, enough for solution of many practical tasks, but it is 
necessary to remark, that the exponential approximation may lead to a significant error in the case of extremely large 
noise intensity, and in the case when the noise intensity is small, the potential is tilted, and the barrier is absent 
(purely dynamical motion slightly modulated by noise perturbations). 
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FIG. 1. Evolution of the nondecay probability for the potential $ (x) = ax 2 — bx 3 for different values of noise intensity; the 
dashed curve denoted as MFPT (mean First Passage Time) represents exponential approximation with MFPT substituted into 
the factor of exponent. 
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FIG. 3. Evolution of the nondecay probability for the potential $(x) = 1 — cos(:r) — ax for different values of noise intensity. 
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FIG. 4. Evolution of the nondecay probability for the potential $(x) = —6a; 3 for different values of noise intensity; the 
dashed curve denoted as MFPT (mean First Passage Time) represents exponential approximation with MFPT substituted into 
the factor of exponent. 
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